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Abstract 

We calculate the probability distribution for the volume of the Universe after slow-roll 
inflation both in the eternal and in the non-eternal regime. Far from the eternal regime the 
probability distribution for the number of e-foldings, defined as one third of the logarithm 
of the volume, is sharply peaked around the number of e-foldings of the classical inflaton 
trajectory. At the transition to the eternal regime this probability is still peaked (with the 
width of order one e-folding) around the average, which gets twice larger at the transition 
point. As one enters the eternal regime the probability for the volume to be finite rapidly 
becomes exponentially small. In addition to developing techniques to study eternal inflation, 
our results allow us to establish the quantum generalization of a recently proposed bound 
on the number of e-foldings in the non-eternal regime: the probability for slow-roll inflation 
to produce a finite volume larger than e^'^^l'^ , where Sds is the de Sitter entropy at the end 
of the inflationary stage, is smaller than the uncertainty due to non-perturbative quantum 
gravity effects. The existence of such a bound provides a consistency check for the idea of 
de Sitter complementarity. 



1 Introduction 



The Universe is accelerating today [1, 2, 3], and it is extremely likely that it was experiencing a 
period of accelerated expansion (inflation) back in the past [3, 4, 5, 6]. In both cases the pressure 
to density ratio is very close to —1 and the local geometry is very close to that of de Sitter (dS) 
space. 

It is plausible that both these periods of inflation are eternal, i.e. some space-time regions keep 
inflating forever. Indeed the most economical explanation for the cosmic acceleration observed 
now is that we are stuck in a metastable vacuum [7], and this results in eternal inflation unless 
the vacuum decay rate F is faster than the expansion rate of the Universe, F > Hq ^. Also there 
is strong evidence that in the past we underwent a phase of inflation driven by a rolling scalar 
fleld, which could have been preceeded by a period of eternal inflation, as predicted in many fleld 
theoretical models [9, 10, 11]. Further, the current picture of the string landscape [12] suggests 
that the observed part of the Universe was created as a result of tunneling from some higher- 
scale eternally- inflating vacuum [13]. All these arguments make the study of eternal inflation very 
important. 

Moreover eternal inflation provides a natural framework to implement Weinberg's solution of 
the cosmological constant problem [14], which is the most plausible so far in spite of the many 
efforts made to flnd an alternative. According to this solution the choice of the vacuum is made, 
at least partially, by anthropic reasons such as the requirement that structures were able to form 
in our Universe. These arguments raise the notoriously difficult and puzzling question of making 
predictions in an eternally inflating Universe. 

On a purely practical side we possess a well developed machinery of quantum fleld theory 
in curved space-time that proved to be very successful in calculating properties of the primordial 
density perturbations with many fine details in the case of non-eternal inflation (see e.g. [15]). The 
applications of these techniques to the case of eternal infiation is instead much more challenging, 
as in this latter regime the size of the quantum fluctuations is large, a non-perturbative treatment 
is required and calculations become in general much more difficult (see e.g. [16] for a recent 
discussion). On the other hand, without a clear understanding of the eternally inflating geometry, 
it might be hopeless to solve the issues raised by eternal inflation such as the measure problem in 
the landscape; this is why we flnd it very important to make explicit and precise calculations in 
this regime. 

On a more theoretical side, dS space appears to share many properties with the black hole 
geometry — most importantly, in both space-times the most natural sets of observers (the asymp- 
totic observers in the black hole case and the comoving ones in the accelerating Universe) see a 
gravitational horizon with the associated thermodynamic properties, such as Hawking temperature 
and Bekenstein entropy, the latter, in the case of de Sitter, being equal to [17] 

where Mpj = 1/Gn with Gn the Newton constant. A flnite entropy suggests that the system is 



^This latter possibility appears rather unlikely and fine-tuned given the non-perturbative nature of the decay, 
although it may receive support from future particle physics data [8] . 
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Figure 1: Set of space-like slices covering both the exterior and the interior of: a) a black hole (left) and 
h) an eternal inflating universe (right). 



described by a finite number of degrees of freedom, or more formally, the Hilbert space describing 
the system at the quantum level has a finite dimensionality equal to e'^'*'^. It is widely believed 
that both black holes and de Sitter space always subsector of a larger theory with an 

infinite dimensional Hilbert space. This is obvious for a black hole in asymptotically Minkowski 
or adS space-times, but less clear for dS space. Indeed, in principle one could imagine a quantum 
gravity theory with, for example, a single positive energy vacuum. However, this is not the case 
in the known string theory landscape [18] and there are general arguments strongly suggesting 
that dS vacua are always metastable with respect to decay to either the Minkowski or the adS 
minima of the potential [19]. This will be the point of view adopted in the current paper and by 
dimensionality of the Hilbert space describing black hole or dS space we always understand the 
dimensionality of the corresponding subsector in a larger, most likely infinite-dimensional, Hilbert 
space (although, it is worth noting that an alternative line of thought is also being pursued [20]). 

Taking seriously the similarity between the causal structures of de Sitter and Schwarzschild 
causes a serious doubt on the validity of the global semiclassical picture of the eternally infiating 
Universe^. Indeed, the remarkable fact about black holes is that the global effective field theory 
description of the space-time claiming to describe both the interior and the exterior of the horizon 
eventually breaks down (see e.g. [23]). 

More concretely, if one considers the set of space-like slices covering both the exterior and the 
interior of a black hole, as shown in fig. la, and insists that the field theory description is valid 
on this set of slices, one comes to the conclusion that the information is lost in the course of the 
black hole evaporation [24]. This conclusion was proven to be wrong [25, 26] by the adS/CFT 
correspondence that provides a description of the system involving black holes (gravity in the bulk 
adS space) in terms of an unitary boundary CFT. Consequently, one is forced to conclude that the 
global effective field theory description breaks down at the time scale of order the evaporation time, 
t ~ R^/Gjy, where Rs is the Schwarzschild radius and Gn the Newton constant (for a recent review 
see for example [22]). After this time the information about the inside observer gets reprocessed 
in the evaporating Hawking quanta, and by insisting on a simultaneous local description of the 
exterior and the interior on longer time-scales one would run into a contradiction with the "no 
quantum xerox" principle (or equivalently with the linearity of quantum mechanics) [27, 28]. This 

^This idea is being pursued by a number of authors, see e.g. [21]. Our discussion mainly follows that of [22]. 
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conclusion is really surprising given that one can always choose a set of slices that avoid the region 
close to the singularity, so that naively one would expect the effective field theory to hold. 

Given the similarity between the causal structures of de Sitter space and that of Schwarzschild 
geometry, one may suggest that also in dS space the global description in terms of a similar set of 
slices (shown in fig. lb) eventually breaks down. Note that these slices are just the conventional 
FRW slices that are commonly used to describe the inflationary Universe. Going on with the 
analogy with the black holes one expects this breakdown to happen at a time-scale of order 
t ~ H~^/Gn or, equivalently, after a period of order S^s e-foldings. One may expect that space- 
time events outside the region containing e'^'*'^ Hubble patches get encoded in de Sitter fluctuations, 
similarly to how the information inside the black hole gets released in the Hawking quanta after 
the evaporation time. 

Note that, unlike in the black hole case, we are not running into any paradox if this does not 
happen, therefore it may well be that pushing the analogy this far is too naive. However, if one 
takes it seriously there is an immediate test to pass. Indeed, in the black hole case, due to the 
presence of a curvature singularity, it is impossible to read the information in the Hawking quanta 
and then jump into the black hole and read the same information again. Similarly, it should not 
be possible to get the same information twice in the de Sitter case as well. 

However, if inflation could last for an arbitrarily long time without becoming eternal, an ob- 
server would be able to flrst read the information in the dS fluctuations and later, as more and 
more space-time becomes visible after inflation, he would be able to access it directly. Conse- 
quently, for the above ideas about de Sitter complementarity to be consistent, there should be a 
limit on how long inflation can last without becoming eternal [22], 

3iV < cSds , (1) 

where N is the number of e-foldings (defined as one third of the logarithm of the total volume 
after inflation) and c is a coefficient of order one. On the other hand, no limit on the number of 
e-foldings is required on those realization that are eternal, because in those cases the observer is 
not able to access all the volume after reheating. 

It was proven in [22] that the bound (1) indeed holds for the classical inflaton trajectory in 
any theory of infiation that does not allow violations of the null energy condition, 

p + p>0, 

where p and p are the energy density and the pressure. On the other hand, violation of the 
bound (1) is possible in theories able to violate the null-energy condition, such as ghost inflation 
[29]. This is actually encouraging and supports arguments establishing the link between horizon 
complementarity and duration of inflation, as also the conventional black hole thermodynamics 
breaks down in theories where the null energy condition can be violated [30, 31]. 

There are two reasons leading to an uncertainty in the numerical value of the coefficient c 
in the bound (1) as proven in [22]. First, at the time when the bound was proposed the exact 
condition for inflation to become eternal was unknown, and it was not even clear whether there is 
a sharp distinction between eternal and non-eternal regimes. This issue was addressed in [16] and 
the conclusion is that there is a sharp transition between these two regimes, with the condition 
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not to have eternal inflation being 

n^^4>l. (2) 

where is the classical velocity of the inflaton field. 

Now it is straightforward to find a bound on the number of e-foldings for the classical inflaton 
trajectory in single-field slow-roll inflation in the non-eternal regime. Namely, one writes 



dN, - Hdt - - ^^^^ i'^) 

where at the last step we made use of the second Friedmann equation. By integrating (3) and 
using the condition (2) for the absence of eternal inflation we obtain (1) with the value c = 1/4, 

3iV, < ^ , (4) 

where Nc is the number of e-foldings on the classical inflaton trajectory. 

However, this does not establish the sharp version of the bound (1) yet. There is another 
reason for the uncertainty in the bound (1). Namely, the above analysis (and that of [22]) is 
restricted to the classical inflaton trajectory. This approximation clearly breaks down when fl 
is of order (but still larger than) one, so that inflation is close to be eternal. In this case, even 
though inflation is not eternal, the typical inflaton trajectory is very different from the classical 
one and can be much longer. More generally, at the quantum level for any value of Q there is 
always a non-vanishing probability for the actual inflaton trajectory to be long enough to violate 
the bound (1) for any value of c. 

In this situation it is natural to study what is the probability distribution for inflaton trajecto- 
ries of different lengths — in other words, what is the probability distribution for the volume of the 
Universe p{V) after inflation. It is not clear a priori what the natural generalization of the bound 
(1) should look like at the quantum level. One might expect that there exists a value of c such 
that the probability to violate (1) is suppressed, for example as non-perturbative quantum gravity 
effects e"'^''^ (which would correspond to p{V) ~ 1/V") or even more, for example exponentially 
with the volume p{V) ~ (which would correspond to order ''^ effects). What we flnd from 
our analysis is that such value of c exists (it is c = 1/2) and that the probability associated to the 
violation of the bound is actually super-exponentially small, i.e. ~ e"'^^'*^. 

To achieve this goal we obtained another result of independent interest. Namely, we calculated 
in an explicit form the probability distribution for the volume of the Universe after slow-roll 
inflation p{V) both in eternal and non-eternal regimes. This offers further insight in the actual 
geometry of the eternally inflating spacetime. While, unlike the density perturbation spectrum, 
this quantity is not of much interest for current observations, it still appears to be one of the natural 
"theoretical" observables to look at in the study of eternally inflating Universes. In particular, 
according to [16], the order parameter for the transition to eternal inflation is the normalization 

/•oo 

Pext = / dVp{V) . (5) 
Jo 



4 



At > 1 this quantity is equal to 1, in agreement with the naive expectation. However, at Q < 1 
the normalization Pcxt becomes smaller then 1, indicating that there is a non-zero probability 
(1 — Pext) for the reheating volume to be infinite, i.e. for infiation to last forever. In this paper 
we will rederive this result in yet another, somewhat more explicit, way. It is also worth noting 
that recently the far exponential tail of the probability p{V) was calculated in [32] in the eternal 
regime. This result was used there to define a "reheating-volume" measure for observables after 
eternal infiation. It appears that the explicit expression for p{V) has good chances to be useful in 
further theoretical studies of de Sitter space and eternal infiation. 

The rest of the paper is organized as follows. We start section 2 with a review of the discrete 
stochastic branching process introduced in [16] to describe infiation. Then we use this model to 
derive a differential equation (32) for the Laplace transform of the probability distribution p{V). 
Similar discrete models were used in [33] and they are essentially equivalent to the stochastic 
description of inflation by Starobinsky [34]. Within the stochastic approach equation (32) is 
known as the "non-linear Fokker-Planck equation" [35] . 

In section 3 we provide approximate solutions for the equation (32) and calculate the proba- 
bility distribution in different regimes. We start with a discussion of the general properties of the 
solutions of (32), and rederive in a new and very explicit way that 1] = 1 is the transition point 
to the eternally inflating regime, i.e that Pext = 1 for ^2 > 1 and Pext < 1 for Q < 1. Then in 
section 3.1 we study the moments of the volume distribution. We show that there is a simple way 
of calculating them without actually solving equation (32) and performing the Laplace transform. 
In agreement with the results of [16] we prove that at Q > 1 sufficiently high moments diverge for 
any value of Q if the inflaton fleld is allowed to take arbitrarily high values. We flnd the values Qn 
such that the n-th moment diverges at fl < fin, and illustrate our method by explicitly calculating 
the average and the variance. 

In section 3.2 we start analyzing the properties of the probability distribution by performing 
the Laplace transform. First, we consider the semiclassical limit ^ 1, and flnd an approximate 
solution for eq. (32) in this limit. It turns out that because of the non-commutativity of the large- 
volume and the large-fl limits this solution does not capture correctly the large volume tail of 
the probabihty distribution where the probability becomes smaller than ~ e~^. The study of this 
solution is also instructive for developing an intuition on how to perform the Laplace transform 
of the solutions of eq. (32). In section 3.3 we apply this intuition for general > 1, when it is 
not possible to flnd an approximate solution to (32) in closed form. By solving this equation in 
different regimes one can obtain enough information to reconstruct the probability distribution in 
the physically relevant case ^ 1. p{V) turns out to be peaked around the average value 

. (6) 

For N < N it takes the following Gaussian form 

p{N) ~ Afe 5^ , n > 1 (7) 

where the width a is equal to 
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(8) 
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Figure 2: Typical shape for the probability distribution of the volume p(y). For small volumes the 
behavior is gaussian with the number of e-foldings (p ~ e~'^^^~^^'^ ); for volumes larger than the average 
value V , piy) follows a power law in the volume (p ^ l/V"-) that eventually turns into an exponential 
law (p ~ e~'^°°^*'^J at large enough volumes (V > Vf,J. When 17 < 1 the exponential tail starts earlier at 
1/ ~ Ve = e'^/(2v^T^). 



As Vt approaches the transition point, f2 = 1, the width a becomes of order one, which is still 
narrow in the regime of large number of e-foldings, 3> 1, the one we are interested in. So 
the most important consequence of the change of Vt is that the average number of e-foldings N 
changes. In agreement with the naive expectation it increases as Vt approaches the transition 
point, but it does not grow a lot: at = 1 the average number of e-foldings is twice as large as 
the classical one. 

At large volumes, N > N, the probability distribution becomes exponential in A^ (or, equiva- 
lently, power-law in the volume V^), 

p{N) oc e-'^^^+v^) = V-^K^+v^) . (9) 

Then we proceed with the eternal inflation regime. First, we study what happens in the 
vicinity of the transition point, i.e. when Q = 1 — e with < e ^ 1. We find that the probability 
distribution p{N) is not changed until A^ ~ vr/ {6^/e). At large volumes, N > n/ (Oy^), it becomes 
much more strongly suppressed, p{N) oc e~'^°'^''*''^ . This behavior is easier to interprete in terms of 
the volume distribution p{V) (rather than the e-folding distribution p{N)). It indicates that if the 
volume gets large enough, V > e'^/*^^^^ the probability for inflation to terminate is exponentially 
small Pext oc e"™"*^*'^. Related to that, we find that the total probability for inflation to terminate 
is of order e~^(^^=) (this also applies for e ~ C(l)), indicating that it is saturated by the small- A^ 
tail of the Gaussian distribution (7). This behavior smoothly matches with yet another regime 
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where one can find the probabihty distribution exphcitly — Q ~ 0. Here the probabihty distribution 
p{V) is exponentially small p{V) oc e~^/^ for all volumes of interest, V ^ 1. 

All these considerations were made in the approximation where the infiaton potential goes 
up to arbitrary high values of the infiaton field. This is clearly unrealistic and we conclude in 
section 3 by discussing what happens in the presence of a "barrier" at large values of the infiaton 
field. As expected, the presence of a barrier affects only the far tail of the volume distribution 
by making it exponentially suppressed at large volumes, p{V) oc e~'^°'^***'^, for any value of fl. 
Note that if the initial value of the infiaton field is not too close to the barrier, this effect is 
relevant only for infiaton trajectories whose probability is smaller than the uncertainty coming 
from non-perturbative quantum gravity effects (~ e"*^"*^). The various behaviors of the probability 
distribution are summarized in fig. 2 that shows the shape of p{V) in the different regimes of the 
volume. 

These results imply that the quantum version of the bound (1) does hold with c = 1/2. In the 
concluding section 4 we give the physical explanation of the behavior of the probability distribution 
and show that an infiaton trajectory with more than S^s/Q e-foldings and such that infiation 
terminates globally in the entire space is super-exponentially improbable. We also speculate on 
the possibility that the value c = 1/2 for the coefficient in (1) that we obtained in our analysis 
might have a natural physical interpretation. In the appendix we cross check our results by 
calculating the average volume directly from the infiaton stochastic equations. 



2 From bacteria to inflation 

As explained in the introduction, we want to calculate and study the probability distribution 
p{V, (f)) of the reheating volume given a certain initial value of the field 0. This calculation does 
not seem to be straightforward, as the only available definition of the distribution is a rather 
formal functional integral formula [16] 



p{V,(f))^ Jv^V[^,(f)\S V-Jd^xe 



(10) 



where X>0 is some vaguely defined measure on the set of all possible space-time realizations of 
the infiaton field, is the probability of a specific realization and tr{x) is the reheating 

time for a given realization as a function of the comoving coordinate x. Evaluating directly the 
functional integral is of course a very hard task. As usual with functional integrals, in order to 
gain more control it is natural to switch to a discretized description of the infiationary dynamics. 
This approach has been recently developed in [16] and similar models have also been studied in 
the context of eternal infiation in [33]. Up to small extensions, section 2.1 is mainly a review of 
the results in [16], which we use to derive the formula for the probability distribution p(y,(f)) in 
section 2.2. The resulting solution for p{V, 4>) will be discussed in section 3. As a cross-check of 
the method, in the appendix we also present a direct computation of the volume average. 
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2.1 Review of bacteria model 



With a biological analogy, consider at t = a bacterium that can live in a discrete set of positions 
along a line (see fig. 3). At t = 1 the bacterium replicates into A^,. copies. Then, each bacterium 
(independently of all the others) hops with probability p to the neighboring site on its right, and 
with probability (1 — p) on the left. Nr and p are fixed numbers. At t = 2 each second-generation 
bacterium reproduces itself, and so on. The analogy with the infiationary system is clear: each 
bacterium represents an Hubble patch; sites are infiaton values. Reproduction is the analogue of 
the Hubble expansion; at every e-folding ~ new Hubble volumes are produced starting from one. 
From then on the infiaton inside each Hubble volume evolves independently, with a combination of 
classical rolling and quantum diffusion. This is represented by the random hopping of our bacteria. 
The difference in the probabilities of moving right and left gives a net drift, and thus corresponds 
to the classical motion. To complete the analogy we have to assume that there is a "reheating" 
site, i = in the figure: when a bacterium ends up there it stops to reproduce and to move 
around — it dies. In the bacteriological analogy the reheating volume corresponds to the number 
of dead bacteria (= non-reproducing Hubble patches) in the asymptotic future. For analogy we 
denote the latter quantity by V, which of course now takes discrete values. Our task is to study 
the probability distribution of as a function of the parameter p. A discrete system like the one 
we described goes under name of branching process, more precisely a multi-type Galton- Watson 
process (see e.g. ref. [36]). 



1-p p 

-o — o — o — — o — o 

i+l i i-1 10 

Figure 3: The branching process. 



To connect to the infiationary case, if is the infiaton, and t is the time of the FRW metric, one 
can make the following identifications in terms of the position j, time-step n, field-space interval 
A0 and time interval At, 

7 = , n = — . 11 

' At ^ ^ 

By taking the continuum limit A(j), At — > in such a way that 

At = ^(A0)^ , (12) 

where H is the Hubble rate in the infiationary process, and by defining as 

Nr = l + SHAt , (13) 
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and p by the relationship 

(l-2p)^ = ^ p=^ + ^/6^^, (14) 

the equation for the probabihty of a bacterium to be at site j at time n becomes the stochastic 
equation for the inflaton [34, 10, 11] 

— a,p(0, t) = -a?p(0, t) + ^^a^p(0, r) , (is) 

where P{(j),t) is the probabihty that the inflaton after a time t has a value (p. In [16] general 
arguments establishing the matching with the continuum limit were presented and checked in 
several calculations performed both in terms of the bacteria model and of the inflaton. 

Let us therefore study in more detail the bacteria model, and consider a branching process 
on a line of length L. A convenient tool to study the branching process is the set of generating 
functions fi"'\sj), where i,j = 0,...,L. These are deflned as power series 



t\s,)^ Y.^,..^A---^^ (16) 

k\...kL 

where pf!^^^ is the probability that, in a branching process that started with a single bacterium 
at the i-th site after n steps, one has /cq bacteria at the zeroth site, ki bacteria at the flrst site, etc. 
It is convenient to combine together all functions /j*-"^ with the same number of steps n into a map 
Fn from the L+ 1-dimensional space of the auxiliary parameters s,j into an L+ 1-dimensional space 
parameterized by the /j's. Also in what follows we will sometimes drop the subscript from the Si 
variables and denote by s a point in the L + 1-dimensional space with coordinates (sq, . . . , s^). 
For example, for a branching process of the sort as described in flg. (3), Fi is given by 

fi,^\s^,...,SL) = So, (17) 
fl\sQ,---,SL) = ((1 -p)s2+PSo)^' , 

/f^(sO, ■■■,Sl) = ((1 -p)Sj+i+pSj_i)^'" , 

/l^(so, = ((1 -p)sL+psL_i)^'" , 

where we have made a speciflc choice of boundary conditions at the site L, which we will refer to 
as the "barrier" . The particular choice of the boundary condition will affect only very marginally 
our results. The main property making generating functions useful is the iterative relation 

= Fi(F„) . (18) 

This property is straightforward to check by making use of the deflnition of the branching process 
and elementary properties of probabilities. 
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Figure 4: Left: Plot of Fi{s) restricted to the hypercube II and with sq = 1, for large p (thick curve). 
The only fixed point in the unit cube is s = 1. Further applications of Fi (thinner curves) drive the curve 
to the Foo = 1 line. Right; For smaller p's a new fixed point Sf enters the unit cube. Now the limiting 
line is F^o = sj. 



We will be interested in the late time behavior of the branching process, which is determined 
by the asymptotic function F^o- The iterative property (18) implies that 

Fi(Foo) = , (19) 

i.e. the set of values of the function F^o is a subset of the fixed points of the function Fi, such that 

F,{s) = s. (20) 

For our purposes it is enough to study the mapping Fi inside the L + 1-dimensional cube Il+i 
of unit size, < Sj < 1. The definition (16) of the generating functions implies that all partial 
derivatives of Fi{s) are positive. Also, the normalization of probabilities implies that 

Fi(i,...,i) = (i,...,i) = r. 

In the bacteria model, it is rather straightforward to see how the transition to the eternal 
inflationary regime happens. For this, it is useful to restrict the function F to the L-dimensional 
hypercube II with sq = 1, which amounts to marginalizing over the number of dead bacteria (see 
eq. (16)). Note that, if the mapping Fi has no other fixed points in the cube II apart from 1 (see 
fig. 4), then 

^ooL,=i = r. (21) 

By definition of the generating functions, eq. (16), this means that in the late time asymptotics 
with probability one there are no alive bacteria at any of the sites. The extinction probability 
is exactly equal to one (inflation ends). The situation changes when a non-trivial flxed point Sf 
solving eq. (20) enters the region 1^ (see fig. 4). Now one has 

^ooL„=i = (l,5/)<l. (22) 
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This implies that, as before, the probabihty to have any finite non-zero number of ahve bacteria 
at any site vanishes. However, the probabihties to have zero bacteria at the various sites, 

pt:U..0-ft\0)^{sf),, i^O, (23) 

are all less than one. This means that there is a non-vanishing probability that the population 
never dies out and that the total number of bacteria grows indefinitely at late times. This corre- 
sponds to the eternal infiation regime. Clearly, this implies that the number of dead bacteria also 
has a finite probability to grow indefinitely; in the context of eternal infiation this translates into 

/■oo 

/ dVp{V,ci>)<l. 
Jo 

In the continuum limit and with infinite barrier this transition happens at Q = 1. 

We would like to stress that F^o is a function of Sq only. Mathematically this follows from 
the convexity of Fi with respect to Sj, i ^ 0, from the fact that Fi is defined between and 1, 
and that Fi{l) = 1 (see fig. 4). This can also be intuitively understood in the following way: 
in the non-eternal regime, in the infinite future there is zero probability of having any bacteria 
alive, so F^ can not depend on any of the Sj, i ^ 0. In the eternal infiation case, the population 
either extinguishes or becomes infinite. In either case, the probability of finding a finite number 
of bacteria at any site is zero, which again forbids any dependence on Sj, i ^ 0. Therefore, 

= Foo(so). Eq. (19) for the fixed point then becomes 

fo°°^ (*o) = So , 

ft\so) = {{l-p)flZ\^o)+pfn\so)Y\ 

ft'Hso) = {{l-p)ft\so)+pftliso)Y^ ■ (24) 

This is the set of equations that determine f^°°\so). Once fi°°\so) is found, one can extract the 
asymptotic probability distribution p^J^^-* using eq. (16). 

2.1.1 An example: the 2-sites case 

We find it instructive to illustrate the general formalism using a simple explicit example. Consider 
the minimal branching process with just two sites and Nr — 2 copies at each reproduction event 
(see fig. 5). In this case the generating functions (16) are particularly simple, 

/d ''('^Oi Sl) = So ; 

fi'\so,s,) = {{l-p)si+psof . (25) 

It is straightforward to apply here the generic results discussed in the last section. fi°°\sQ) is 
given by the fixed point of the mapping defined in eq. (25) inside the unit interval < Sj < 1. We 
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Figure 5: The 2-site branching process. 



obtain 



flr\so) = so. (26) 
.(oo). 1 - 2pso + 2p2so - ^1 - 4p(l - p)sq 

fi (^o)- . (27) 

The extinction probability is given by 

P -V„ _ .(oo)..._ l-2p + 2p^-v/(r^2^ _f 1 p>l/2 

where pk = Pi^\- The extinction probabihty indeed drops below one for p < 1/2. We get the 

quite intuitive result that the critical probability is pc = 1/2. By Taylor expanding f^'^\so), we 
can obtain the probability of the volume pk as in eq. (16) 



which allows us to extract the large- /c asymptotic of the probability distribution pk, 

Pk ^ Lg-'^^^g^W^rt) . (30) 

We see that the probability for large volumes goes to zero exponentially fast for any p ^ pc = 1/2 
(notice that Ap{l — p) < I for any p, and becomes equal to one only for p = 1/2). At p = pc the 
exponential suppression disappears and we are left with a power law behavior, corresponding to the 
singularity of /i (sq) at Sq = 1, as a result the probability distribution becomes not normalized 
for p < Pc- This can be understood in the following way. So far we have studied the probability 
distribution in the infinite time limit. If instead we let the process go on only for a finite number 
of time steps n, the probability distribution develops a bump at large volume (see [16] for details), 
which can be roughly thought of as a piece proportional to 0(1 — Pcxt)S{k — 2"). For p < pc, in 
the limit n — > oo, the 5— function migrates to infinity and makes all the moments of the volume 
diverge. In our calculation, we can not see such a 5-function because we have already taken the 
limit of infinite time steps, but still we are left with a non-normalized distribution. Because of 
this even for p < p^ all the moments of the volume diverge even if the probability distribution 
goes exponentially fast to zero at infinity,. 
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2.2 The equation for p{V) 

We want now to take the continuum hmit of the bacteria model as explained in the beginning of 
section 2.1. According to eq. (11), we define f^°°\(l); So) — f^A4>i^o) so that the system in eq. (24) 
now becomes 

/(-)(0,so) = So, (31) 
/(-)(0;so) = {{l-p)f^^\<f> + A<l>;so)+pf^\<t>-A<t>;so)f\ 



f^\<Pb;so) = ((l-j9)/(°^)(0,;so)+p/(°°H0.-A0;5o)) 



Nr 



where 4>h represents the value of 4> at the barrier. 

By taking the limit A0 — > At — > according to the prescription given in eqs. (12), (13), and 
(14) we obtain the following second order differential equation 

^o) - ^^/^-n^; ^o) + ^o) log [f-\<l>; .0)] = , (32) 

with the following two boundary conditions 

/(°^)(0;so) = so, (33) 
= 0. 
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in agreement with [35, 32] where the same equation has been obtained in a different way. The 
second derivative term in eq. (32) comes from the quantum fluctuations of the inflaton (equivalent 
to the random jumps of the bacteria), the first derivative term comes from the classical drift (in 
fact it is proportional to Q), and the log in the last term comes from the production of independent 
Hubble patches (equivalent to the reproduction of bacteria). 

This equation will be central in the rest of paper: /(0; sq) is the Laplace transform of the 
probability of obtaining a certain reheating volume starting from any initial value of 0. Indeed, 
in the discrete model, from the definition in eq. (16), the generating function at infinite time 
is connected to the probability pj^k of having k dead bacteria at infinite time starting with one 
bacterium at the site j by 

oo 

ft\'o)^J2PJ^ksl (34) 

k=0 

Taking the continuum limit we get 

POO 

f^\'P;so)= / dVp{<l>,V)s^. (35) 
This expression can be inverted to obtain the probability distribution for the volume 

1 /■7+ioo 

Pi<P,V) = — / d(-log(5o))/(°^)(0;5o)e-^^°^(^°) (36) 
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where 7 must be chosen such that Re(7) > — log(sQ"^^) for any singularity Sq°^ of sq) ^. Notice 
that in eq. (36) we need to analytically continue the function f^°°\(p; sq) to unphysical values of 
-So ("So ^ [0, 1]), which actually dominate the integral at large volumes as we will see in section 3. 
We have therefore obtained a procedure to compute the probability distribution of the reheating 
volume: solve the differential equation (32), and then perform the anti-Laplace transform (36). 

3 Probability distribution of the volume after inflation 

In the previous section we saw that the probability distribution of the volume can be calculated 
in two steps. The first is to solve the differential equation (32). For convenience we rewrite it here 
as 

fir; z) - 2v^/(t; z) + fir; z) log[/(T; ^)] = , (37) 

where the dot represents a partial derivative with respect to r and f{T;z), r and z are related to 
j{°o)(0; So), (f) and sq of the previous section via 

f{r;z) ^ /(°°)(0;so), (38) 
T = 2nV6^ = eVnNc , (39) 
z = -log(so), (40) 

with Nc being the classical number of e-foldings {Nc = —H(p/(l)). The solution /(r, z) has also to 
satisfy the following boundary conditions 

/(0;^) = so = e-^ (41) 

/(r6;^) = 0, (42) 

and the constraint 

/(r; z) e [0, 1] . (43) 
The second step is to calculate the integral 

r>0++joo 



P(^,^) = ^/ dz, 
2m Jo+-ioo 



f{r;z)e^\ (44) 



that gives the probability distribution p{V, r) to find a volume V at the end of an inflationary 
phase that started with the infiaton at the position = i?r/(27r\/6). Recall that the volume 
V that enters in eq. (44) is dimensionless because it has been rescaled by the initial volume Vq 
before inflation, i.e. V = VoI/Vq, with Vol being the physical volume. Notice that some properties 
of p{V) can be obtained without evaluating the integral (44), but just by studying the solution 



^From eq. (35) we see that /^°° ■*((/>; So) cannot have singularities for Re(log(so)) < 0, therefore eq. (36) holds for 
every 7 > 0. 
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Figure 6: The potential U{f) of eq. (47)- limit ti, oo, the solution of the mechanical problem 

(37) represents the motion of a particle starting at t = in f = sq, rolling uphill with an anti-friction 
term and stopping on the top f = in an infinite time. 



f{T; z) around the point z = {]. Indeed the momenta of the distribution are simply related to the 
Taylor coefficients of /(r; z) around z = since 

(0 = dvvy{y,r) = {-ir 

It follows that the total probability to exit inflation globally is just fixed by /(r; 0), 

/■oo 

P,^,= / dVp{y,r) = f{T-(i). (46) 

We proceed now with the study of the solutions of eq. (37). This equation describes the motion 
of the particle in a potential 

W) = ^(logf -1) , (47) 

with an anti- friction term (see fig. 6). Unfortunately an explicit solution to eq. (37) is not known. 
However we will still be able to recover many properties of p(V, r) by analyzing the analytic 
structure of the solution and by studying the problem in different limits. The boundary conditions 
(41) and (42) as well as the condition / G [0, 1] constrain the solution to start from the point 
/ = So at r = and to travel up-hill up to some point /=„ G (0, sq) where the velocity reaches zero 
at time r^. We will focus on the case where +oo, and come back to the case of finite at 
the end of this section. In the r;, +oo case, the solution has to reach zero velocity only after an 
infinite time, and in doing so, it has to stay always within the interval / G [0, 1]. The only way to 
achieve this is for the solution to reach the top of the hill / = in infinite time. 



(45) 

2=n 
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Let us see that such a solution exists for any sq < 1. Clearly for large enough initial velocities 
the sohition overshoots the maximum and goes to the / < region, while for small enough veloc- 
ities, the solution does not reach the top. Therefore, there is a critical initial velocity separating 
these two regimes such that the solution stops at / = in an infinite time. This can also be 
explicitly verified by finding, as we will show later, an asymptotic form of this solution in the 
region < / ^ 1, 

/ ^ e-2("+"°)' , (48) 

where Tq is an integration constant. 

The case sq = 1 is different. For a solution that starts at / = 1 at early times we can 
approximate the potential U (/) with an harmonic oscillator to obtain 

/-2v^/ + /-l = 0, (49) 

whose general solution is of the form 

/ = 1 - (Ae^^^ + Se-^^^) . (50) 

This shows that for ^2 > 1 the solution is 'ovcr-anti-dumped' and docs not have a turning point, 
while for ^2 < 1 the solution can oscillate. This behavior persists at the non-linear level as well — it 
is straightforward to check that the turning force due to the potential (47) is always smaller than 
the turning force due to the harmonic potential in (49), so that no turning point exists for fl > 1 
for the non-linear mechanical problem (37) as well. Consequently, at Sq = 1 the solution that 
stops on the top of the hill exists only at Q < 1. This solution describes a non-trivial fixed point 
(22). Its presence indicates that inflation is eternal at < 1. Instead, at > 1 the solution 
that reaches the top of the hill in an infinite time r — > -|-oo starts at / = 1 in the infinite past, 
T = —oo; while all solutions that start at / = 1 at finite time overshoot the top of the hill. 

We illustrate the behavior of the solutions in the two different regimes in fig. 7. This plot 
makes very explicit the transition to the eternal regime aX Q = 1. At ^2 > 1 by taking the limit 
So — 1 one gets /(r; 0) = 1 for every r so that the extinction probability is Pext = 1- On the 
other hand, in the same limit at Q < 1 one obtains a non-trivial function /(r; 0) that leaves the 
origin in a time of order 1 / Vl — ^, therefore there is a non- vanishing probability to inflate forever, 
nxt = l-/(r;0). 

In the following subsections we will present estimates for p{V, r) in different regimes and discuss 
finite-barrier effects. However, before entering the analysis of the probability distribution itself, 
we would like to discuss first its moments, which can be derived exactly in a straightforward way 
from the master differential equation (37). 

3.1 Moments of the probability distribution and critical points 

In spite of the fact that the phase transition happens at the critical value Q = 1, in [16] it was 
shown that, in the infinite barrier limit, the moments of the probability distribution start diverging 
at different values of Q: the higher the moment, the higher the value of Q at which they diverge. 
On the other hand, in the finite barrier case all the moments diverge at the critical value ^2 = 1. 
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Figure 7: Schematic plots of the solutions f{T;z) as a function of t for different choices of the boundary 
condition z and for > 1 (top) and 0, < 1 (bottom). For ^1 > 1 when z ^ the solution approaches 
/(t;0) = 1 for every finite t's — the probability to end inflation globally is 1; for i7 < 1 m the same 
limit the solution approaches a non-trivial function with f{T;0) ^ 1 for large r — the probability to end 
inflation globally is small. 



In the first part of this section, we will rederive these results in a rather quick and different way, 
obtaining also the general formula for the critical value Q at which each moment diverge. Then, 
in the second part, we will describe a procedure derived from the master differential equation (37) 
that allows us to easily compute the explicit value of each moment as a function of fl. 

As shown in the previous section from eq. (45), it is possible to extract the moments of the 
probability distribution p{V, r) directly from /(r; z) without the need of performing explicitly the 
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anti-Laplacc transform. The n-th moment is indeed proportional to the n-th derivative of /(r; z) 
with respect to z, at z = 0. Divergences in the moments thus correspond to a non-analyticity of 
/(r; z) at z = 0. Therefore it is enough to study the solution of the differential equation near the 
point z — 0. Notice that, before entering the eternal regime, at Q > 1, for every finite fixed value 
of r and for smaller and smaller values of z, the solution is better and better described by the 
linear approximation of eq. (50) (see fig. 7), which we conveniently rewrite here as 

/iin(T; z)^l- e'"-(^+^°) - (7e'"+(^+^°) , (51) 

where 

and a and tq are the two constant of integration. In the infinite barrier case the boundary condition 
(42) corresponds to requiring that the solution stops on top of the hill (/ = 0) in an infinite time. 

Notice that since both this boundary condition and the differential equation are invariant under 
shifts of T we are left with a one-parameter family of solutions that are related by a shift in r, i.e. 
the parameter Tq. The latter is fixed by imposing the boundary condition (41), namely 

/ii„(0; z)^l- 6""-^° - (7e^+^° = e-' . (52) 

This means that all the dependence on z is in the parameter tq, while a is fixed by the boundary 
condition at infinity and is independent of z. Using eq. (52) the solution can be rewritten as 

/ii^(r; z) = 1 - (1 - e-') e'""^ - ae'^+^'> (e'"+^ - e'""^) 

~ 1 - ze""-^ - az<e'^+^ , (53) 



where in the second line we have used the approximate solution tq ~ log(^)/a;_ from eq. (52) and 
dropped a subleading exponent in the last term. 

Notice also that the last term, in general, is not analytic in 2; = 0. Indeed if we calculate the 
n-th. derivative of /iin(T; z) with respect to z we get 

fi::!{r;z)r^z-l-e--^ (54) 
so that the moment (l^") starts diverging when uj\ becomes smaller than n, i.e. at 

Q^(!i±i£. (55) 
An 

For instance, in the case of the variance (n = 2), we get Vt = 9/8 as critical value, in perfect 
agreement with the lengthy calculation of [16] that used directly the infiaton stochastic equations. 
In the presence of a finite barrier this argument breaks down because now also a depends on z 
and the analytic structure around the origin changes. In fact, for small enough values of z, f{T; z) 
is still well described by the linear approximation around f — 1 even at the barrier t;,. It is not 
difficult then to solve the differential equation in the linearized limit with a finite barrier. What 
we get in this case is 

hn{T;z) = 1- {1-e-') ± , (56) 
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which is analytic in z = for all values of fi, implying that all the moments converge for 1] > 1. 
Although we expect the linearized approximation to work better and better as z ^ 0, we cannot 
be sure yet whether non-analytic terms may arise from subleading non-linear corrections. Still this 
argument suggests that the finite barrier case behaves differently than the infinite barrier case. 

As a proof of this statement we will now present a method to derive exactly all the moments. 
Indeed, even though we are not able to solve analytically the non-linear differential equation (37), 
the equations for the moments are linear and can be solved explicitly. They can be obtained by 
simply deriving n times the eq. (37) with respect to z at z = 0. For example by deriving (37) 
once with respect to z we get 

f_2Vnf + f + flogf = 0, (57) 

where "dots" represent derivatives with respect to r and " ' " with respect to z. Since for z — 
/ = 1 we get a linear differential equation for /q = /'(r; 0) = —{V) with solution 

= Ae^+^ + Be''-^ . (58) 

The constants of integration A and B can be fixed using the derivative of the boundary conditions 
(41) and (42), namely 

/^(0) = -l, /^(t,) = 0. (59) 

This way we get 

{V) = -/^(r) = ^ , (60) 

which in the large limit gives (see fig. (8)) 



lim (y) = e'^-^ = e '^+V^^ . (61) 

This results nicely agrees with the explicit calculation from the inflaton equation (see the ap- 
pendix) and with the result from the probability distribution p{V, r) that we will derive in the 
next sections. Notice that for large Q one recovers the classical limit for the average volume 
Vc — e^^=. With very httle effort eq. (37) gave us the formula for the average volume in both finite 
and infinite barrier cases. 

Roughly with the same amount of work we can obtain also the expression for any higher 
moment. By deriving eq. (57) with respect to 2; a second time we obtain 

f" - 2y/Uf" + /" + /" log / + ^ = , (62) 

which again gives a linear differential equation at 2; = 0, 

- 2v^/^ + = -f^' , (63) 

this time with a non-homogeneous source term. The latter however is an exponential, therefore 
the differential equation can be easily solved analytically. Imposing the boundary conditions 

/^'(0) = 1, f^{n)^o, (64) 
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Figure 8: Average volume as a function of Q, with = ^ff^ = 



the result for the second moment is 



+ 



+ ■ : : ^ . . - + 



^+2™^ 



2 ' (65) 

{2uj\ - 1) (e^''/'^+ - e^b'^+cj^) 

where the length of the expression indicates how non-trivial it would have been to obtain this result 
directly from the inflaton equation. In the large barrier limit the asymptotic form of eq. (65) for 
r2 > 1 reads 

/T/2\ ^+ ( 1 _ 2^^^ e^^-^ + 8K - + 1) -K~2V_r,+c.4-r /aa^ 

^ ' col J ^i(2a;+-l)(2-^2)e , l^oj 

the divergence at tu^ = 2 (i.e. Q = 9/8) is manifest, in particular the last term vanishes for 
a;^ > 2 in the limit Tb oo, while it explodes for < 2. It is less manifest from eq. (65) the 
fact that, with a finite barrier, there is no divergence; for this propose we give here the expression 
for eq. (65) at = 9/8 

(^')L=9/8 = - e-^/^)e^^ + . . . , (67) 
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which shows that, up to sub-dominant terms in the large barrier hmit (the dots), the result is 
finite but linear in Tf,, explaining the divergence in the infinite barrier case. 

We could keep going calculating higher moments, indeed for the n-th moment we have just to 
solve the following linear differential equation 

ft^ - 2VUft^ + = J(") , (68) 

/o^"^(0) = (-!)", (69) 
/o^'^^(r.) = 0, (70) 
where the source J^"^ is a polynomium of degree n of the lower moments {f^''^ with k <n) 

j(") = 9r[/(log/-l)]L^o, (71) 

which will then be a sum of exponentials of the form e^^^'^ at most of degree k = n. 

Iterating the analysis it is possible to check that, in the infinite barrier limit, the critical value 
of Q where the n-th moment diverges perfectly agrees with eq. (55), while for finite barriers the 
moments converge for every Q > 1. 

3.2 > 1: Classical limit 

Let us now find approximations for the probability distribution p{V)m different regimes by directly 
using the Laplace transform formula (44). The first case we will study is the limit Q ^ 1, which 
corresponds to the conventional slow-roll inflation, far from the eternal regime. As expected, we 
will see explicitly that in this case the volume probability is sharply peaked around the volume 
corresponding to the classical inflaton trajectory. The main reason to start with this case is that 
we will be able to find an explicit solution to the mechanical problem (37). This will help us 
to develop an intuition on how to proceed also for generic values of ft, where such a solution is 
absent. 

To analyze the large fl limit it is convenient to rescale the time variable r as 

T = 2Vnf . (72) 
The new time variable f measures the number of classical e-foldings N^, 

f = 37Vc- (73) 
Then the mechanical equation (37) takes the following form 

It is useful to present / in the exponential form 

/ = e-^ (75) 
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then the function g satisfies 



1 

4Q 



df^ \ df 



1 + 5 = 0. (76) 



The most straightforward way to do the large Q expansion would be, as a zeroth order approxi- 
mation, to drop the first two terms in this equation. This would give 

g = e^+^° , (77) 

where fo is an integration constant. This translates into 

/ = e— % (78) 

when we fix tq to match the boundary conditions for /. The corresponding probability distribution 
that we obtain from eq. (44) is 

p(y,r)=5(y-e^^=), (79) 

i.e. the infiaton follows the classical trajectory, exactly what expected in the classical limit where 
quantum flTictuations can be neglected. However, there is a problem to use this solution as a basis 
for a systematic expansion around 1/Q = because the corrections due to the dropped terms in 
(76) are not always small. Indeed, for f 3> logfl one has g^/Q ^ g,g in this case. We can get 
around this problem by keeping this dangerous term in (76), and dropping only the very first one 
(similarly to the WKB approximation). Now we have the following equation 

that after integration gives as a solution of the following algebraic equation, 

Ge^ = e^+^° , (81) 

where 

For g/fl ^ 1 this gives back the previous result (77), while in the opposite limit g/fl ^ 1 one 
gets 

g^n(f + fo)^. (82) 

Importantly, for this solution the dropped term g/ft is small compared to the other ones in (76) 
in both limits as long as n|G + 1| 3> 1. Therefore this solution provides a basis for a consistent 
1/Q expansion which is valid everywhere apart in the small region IG + 1| < 

Let us see what probability distribution one gets from the above solution. We can find fo{z) 
by imposing the boundary condition 

/(O; z) = e-' = e-^(^=°) , (83) 
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which gives 

where we have defined 

Then expression (44) for the probabihty distribution takes the following form 

fO++ioo 



l-Kl Jo+-ioo 



(86) 



where G is a function of x via eqs. (81), (84) and (85). We can try to evaluate this integral using 
the saddle point approximation (we will check later whether this is a good approximation). The 
result reads 

p{V, r) ^ , ^ e-^("°) , (87) 

where S{x) is given by 

S{x) = Q (G^ + 2(7 - xV) , (88) 
and Xq is the saddle point satisfying the equation 



S'{xo)^n(^2{l + G)G'^-V^ 



= . (89) 



X=Xo 



By taking the derivative of (81), and plugging in the resulting expression for G' into the saddle 
point condition (89) one gets 

G = VA (90) 
at the saddle point. If we substitute this expression back in eq. (81) we get 



Ge^ = Al^e^*^ =>le"+^, (91) 
1 /e^^-\ 



V-1 \ V ^ 

Note that for V < Vc = e^^" the value of A is positive and the saddle point is at real and positive 
z (see fig. 9). On the other hand for V > Vc the value of A is negative and the saddle point 
is located at negative z To understand when the saddle point approximation is applicable for 
calculating the integral (86) it is important to note that our approximate solution has a branch 
cut starting at the point Zcut ~ — 2Qe~'^~^, where G = — 1 (equivalently, g — — Q). An easy way 
to see this is to use eq. (80) as an equation for dg/df, namely 



df V V n 



^'=20f-l + ,/l^) . (93) 



■^This implies that for V > Vc the saddle point corresponds to a value for the boundary condition z that is 
outside of the physical region z > 0. We are using the analytic continuation of the solution, as it is implicit in the 
definition of the inverse-Laplace transform of eq. (86). 
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Figure 9: Contour of the anti-Laplace transform integral in the z-plane (on the left) and probability 
distribution of the volume (on the right) for O 3> 1. For small volumes the integral can be solved via 
saddle-point approximation and gives a gaussian-like distribution (in blue). Near V = eVc the saddle point 
hits the cut, the contour of integration can be deformed around the cut and for V > eVc the distribution 
becomes exponentially small (in red). The large limit does not capture the cut between z^ut o,nd z = 0, 
which deforms the tail of the distribution (from a value of V ^ {Vc,eVc) where p{V) ~ e~^) making it 
power-like in the volume (in yellow). 



We see that a.t g = —Q there is a discontinuity in the imaginary part of dg / df along the real axis 
in the ^^-plane at g < —Q, leading to a cut for g a.s a function of z at real z < z^ut- 

From (90) and (92) we find that the saddle point hits the cut at ~ eVc. As long as the 
saddle point does not hit the cut, i.e. for V < e\4, we can use the saddle point approximation to 
perform the integral. Plugging the solution in eqs. (90) and (92) back into S we get 



S{xo) 
S"{xo) 



V 



V 



V-1 



log 



log 



v-1 



log 



, Vc 



V-l 



and thus 



where the prefactor J\f is equal to 



V<eV,, 



1/2 



(94) 
(95) 

(96) 
(97) 



Notice that ^"(xo) is large for exponentially large volumes, making our saddle point approximation 
justified. 
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For large volumes {V ^ 1) the formula (96) simplifies to 



le-^(i°^^)\ V<eV,, (98) 
which we can also rewrite as the probability distribution to have N e-foldings 

m Nc) = SVpiV, r) 3^ ^|3iV - i j ^-nisN-sN^f ^ < _ (gg) 

This distribution is a gaussian centered around the classical number of e-foldings iVg (see fig. 9). 
The spread is of order 1 / and goes to zero as O goes to infinity reproducing the 5-function of 
the classical limit. 

When the volume becomes approximately eVc the saddle point reaches the cut and we cannot 
perform the saddle point approximation anymore. In this regime, we can still close the contour of 
integration on the left around the cut (see fig. 9), and perform the integral along the discontinuity. 
We obtain 

1 r + oo 

P{V, r)^— d\z\ 2i Im[/(f; -|^|)]e-^l^l . (100) 

It is straightforward to verify the |/(f; z)\ does not grow faster than el^l at large z. Therefore, at 
large volumes the integral above is dominated by values of z very close to the cut, with a spread 
of the order 1/V. We can thus write approximately 

p{V, f) ^ e^-*^ ~ 6-^9^/^= , V > eV, , (101) 

where we have ignored power corrections in the volume. We see that the probability distribution 
has an exponential tail in V that starts many standard deviations away from the average Vc- This 
result confirms that after the saddle point hits the cut one cannot use it any longer to evaluate 
the Laplace transform. Indeed, if one keeps using the saddle point approximation blindly one 
would obtain the gaiissian behavior for the probability distribution up to arbitrary large volumes, 
in contradiction to (101). 

However, there is a problem with the behavior (101) as well. Namely, this result disagrees with 
our analysis in section 3.1, where we found that high enough moments of the volume distribution 
diverge at any value of fi. Related to this, we proved there that the function / has a cut starting 
at ^ = 0, while here we are finding the origin of the cut at z = Zcut < 0. 

The origin of this discrepancy is the non-commutativity of the large Q limit and the large 
volume limit. One indication of the problem is that the large Q expansion breaks in the vicinity 
of G = — 1 — precisely where the cut of the approximate solution starts. Even more relevant is the 
following observation. The leading non-analytic term in (53), that gives rise to the cut starting 
at z — 0, is proportional to z*^. At small z this term is non-perturbatively small in the large 
fl expansion. However, it dominates the behavior of p{V) at large volumes and gives rise to the 
power- law tail proportional to V'"^^. Our approximate solution misses the corresponding part of 
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the cut between and Zcut- Note, however, that this part of the cut becomes important only at 
volumes much larger than the average, where the probability is exponentially suppressed in Q, 
p{V) oc e^^. Consequently, our approximate solution correctly reproduces the shape of p{V) up 
toV< eVc- At larger volumes instead of the exponential behavior (101) one gets the power-law 
tail. In the next section we will discuss this tail in more details for general Q > 1. 

3.3 Q > 1: Approaching the phase transition 

Let us now reconstruct the probability distribution of the volume for generic > 1. Unlike in 
the previous case, we do not have any small parameter that allows us to find an approximate full 
solution to the mechanical problem (37). However we will be able to find an approximate form of 
the probability distribution p{V) practically at all V. 

In order to do this, we notice that we can solve the differential equation (37) (and equivalently 
eq. (74)) in two different regimes. When / ~ 1 (equivalent to <^ 1), we can linearize the 
potential and obtain the solution 

/iin(T; z) = l- e"-(^+^) - (7e'^+(^+^°) , (102) 

where u!± = VTl, ± — 1 and a and Tq are integration constants. For the linear approximation 
to hold, it is enough to impose that r + tq ^ —1. Then, independently of the value of a, for large 
enough r + tq we can approximate the solution as 

/lin ~ 1 - e'^-(^+^°) . (103) 

The second regime in which we can solve the differential equation is when / — (equivalent to 
g ^ 1). In this regime the two dominant terms in (76) are those proportional to g and (|^)^; by 
dropping the other terms one obtains 

(t + ti )2 

f^h- e ^ , (104) 

where Ti is an integration constant. By plugging this solution back into the equation, one may 
check that this approximation indeed holds as long as g ^ 1, i.e. for |r + ri| ^ 1. 

Notice that the constants of integrations in both cases can be absorbed into a shift (tq or Ti) 
of the "time" variable r. Tq and Ti are in general not equal — they differ by an unknown constant 
of order one set by the matching of the two solutions. However, later we will be interested in the 
large (r + Tq) limit where such a constant can be neglected and the Tq and Ti can be taken as 
equal. 

Let us sec now that this information is enough to reconstruct p{V) almost for all V. First of 
all, we need to know how / depends on z. From section 3.1 (see eq. (53)), we know that /(r; z) 
has a branch cut a.t z = 0. We could perform the integral by closing the contour around the cut, 
or by using the saddle point approximation. Let us start by seeing if and where we can use the 
saddle point approximation. 

For large volumes the saddle point is expected to lie at small 2; (/ ^ 1), where the exponential 
suppression in (44) is milder. In this case, near r = 0, the linearized approximation holds, and we 
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can use /un to relate tq to z. Next, in order to evaluate the integral (44), we need to impose that 
the saddle point belongs to one of the regions where our asymptotic solution works. Within the 
linearized regime there is no saddle point, so we will check whether the saddle point exists in the 
gaussian region at r + Tq ^ 1. 

Let us see how far the oTitlined procedure takes us and let us begin to implement it. Assuming 
that z at the saddle point is small (we will check this assumption later), we can determine Tq from 
the boundary condition 

e-' ^ fUr ^ 0; z) ^ 1 - e'^-^^ To^—logz. (105) 



00- 



We can now substitute this value for tq in fg and perform the integral (44) obtaining 

1 



and S{z) is given by 
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We see that for relatively small volumes Zq is positive and small, which justifies our assumption 
to use the linearized limit to match tq with z. As V grows, zq moves towards zero (see fig. 10) 
and reaches the region where the gaussian approximation breaks when V V, with V given by 
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e — e 
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(109) 



Even though we can not trust the gaussian approximation for y ~ y, we can try to follow the 
location of the saddle point, and we can see that it moves to negative values for V >V, reaching 
the location of the cut. This further justifies the approach we will take in the regime V >V, that 
is to do the integral along the cut. 

For the moment instead let us concentrate on the regime of volumes V ^ V, where we can 
apply the saddle point approximation. Substituting the value of the saddle point back into S{zq) 
we get 



S{zo)^^(T-^\ogV 



3Nr - 3N 



^ 2 



(110) 



We see that r + tq at the saddle point is large whenever S{zo) is large, i.e. at large and Nc as 
long as A^ < 2(f2 — ^^2(^2 — l))A^c, (compatibly with the condition ^ ^), so that in the same 
regions we can trust the use of fg for the saddle point. The corresponding value of S"{xo) is 
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Figure 10: Contour of the anti-Laplace transform integral in the z-plane (on the left) and probability 
distribution of the volume (on the right) for $7 > 1. For small volumes the integral can be solved via 
saddle-point approximation and gives a gaussian-like distribution (in blue). Around V = V (gray region) 
the saddle point hit the cut, which starts at z = 0, and our approximations for the solution of the 
differential equation break down. At larger volumes there is no saddle point anymore but the contour 
of integration can be deformed around the cut, giving a distribution tail that follows a power-law in the 
volume (in red). 



As in the large Q limit, S"{xo) is large for large volumes. So we conclude that the probability 
distributions for the volume after inflation for V and for generic > 1 has the form 

p(\/,r) ^Are-^K^+v^)'M^)]' =Are-^[¥(i+v^)-3A^^]' , y <v , (112) 

where in the last step we used that 3A^ = log V. We can trust this expression as long as log V > 
N > 1 (see fig. 10). _ 

Let us now deal with the regime V >V. In this case, we have seen above that the saddle point 
Zq enters the region very close to zero where we can not trust anymore the gaussian approximation. 
In section 3.2, eq. (53), we saw that /(r; z) has a branch cut at the point z = 0. Near the branch 
cut, for sufficiently small values of z, we can take the linear approximation /iin(T; z), which reads 

fnn{r; z) = 1 - (1 - e") e"""" - (Te'^+"« (6"^+" - e"^"") 

~ 1 - ze'^-^ - (rz^e'^+^ , (113) 

where in the last line we have used the approximate solution tq ^ log(2;)/ci;_ of eq. (105). It is 
important to keep in mind that a does not depend on z, but it is only fixed by the boundary 
condition at r ^ +oo. 

Since uj^/u^ = oj\ in general is not integer, from eq. (113) we see that /iin(T; z^ has a branch 
cut at z = 0. In the ^ 1 case we did not see this cut starting at z = because, at large fi. 
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ujj^/oo^ ~ 4^2. At small 2;, this is a singularity that appears only non-perturbatively in 1/0 and 
could not be seen in a perturbative expansion in l/Vt as we did in the previous section. 

We can now perform the integral of the discontinuity along the cut. As in section 3.2, we have 
to compute the imaginary part of /(r; z) along the cut and then integrate it 

1 r+co 

piV, r)^^.J^ d\z\ 2i Im[/(r; -k|)]e-^l^l . (114) 

Since /(r; z)e^^ rapidly decreases when the real part of z is negative, at large enough volumes 
there is an interesting regime where the integral is dominated by small values of z, such that we 
can perform the integral using the linearized expression for /(r; z). In this regime we have 

Im[/u„(T;^)]eut~e'^+^z^ . (115) 

where we neglected order one coefficients. By performing the integral (114) we obtain 

f± 

P(^,r)~i(^y" for ^>F^, (116) 

where we used that e'^^^ = V . The condition on the right, that determines how large the volume 
should be for this approximation to work, comes from imposing the validity of the linear approx- 
imation for f{T;z). Indeed, the integral is dominated by values of z around Zs — For the 
linearized approximation to work funizsir) should be close to one. By plugging in the value for 
Zs in the expression for /un, we find 

+ , (117) 



which implies that the integral on the cut is well approximated by the integral of the linear solution 
for V > ^V. 

Similarly to the classical limit, the distribution for V ^ V is a, gaussian in the number of 
e-folding centered at V with still a quite narrow width, of order one e-folding. Our method does 
not allow to reconstruct p{V) in the vicinity of the average, V < V < Va;+/a;_. Note, however, 
that close to the eternal regime is of order one, so that this range of volumes is not big. 

On the other hand at large fl, when uj^/uj- is also large, we were able to find the probability 
distribution up to K ^ eV, where p{V) was already exponentially suppressed as e~^. 

For V > Vuj^/uj- the tail of the distribution in V is power law. 



p{v,t) ~ y"^"^ ~ 



This tail agrees with our results in section 3.1 on the divergence of the multipoles (55). 

As a cross-check of these results note that they imply that the average volume after inflation 
is given by 



{V) c^e'+V'-h , (118) 

in agreement with the direct computations in section 3.1 (eq. (61)) and in the appendix (eq. (175)). 
We see that as fl approaches the critical point the average number of e-foldings shifts from Nc in 



the limit Q ^ 00 to 2A^c in the fl ^ 1 limit. 
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3.4 < 1: Inside eternal inflation 

We now begin to explore the regime of eternal inflation. As Q, crosses 1, the solution of eq. (37) 
changes its form. The behavior of /(r; z) around / ~ 1 is not overdumped anymore. As discussed 
before, the normalization of the probability distribution drops below 1 in this regime. 

We would like to follow this transition carefully. To this purpose, we take Q = 1 — e with 
< e <^ 1 and study the volume probability distribution within the eternal inflation regime by 
expanding in e. We will follow the same strategy as in the previous two subsections. We expect 
that the solution /(r; z) has a branch cut in the complex z plane, which allows us to perform the 
inverse-Laplace transform either with a saddle point approximation, if the saddle point is away 
from the cut, or along the cut itself. 

If we decide to apply the saddle point approximation, we can concentrate on large enough r, 
so that the saddle point lies in the region / — > where the solution is well approximated by 

fg{T;z)^e T + To»l, (119) 

(we will check later what is the corresponding range of volumes) and tq can be determined in 
terms of z in the linearized regime (as long as 2; <S 1, as we will check below). The linearized 
solution is now given by 

/iin(r; z)^l- (7e^(^+^°) cos (Vn^ir + tq)) « 1 - ae^+^° cos (v^(r + tq)) . (120) 

Note the oscillatory behavior of this linearized solution. By rescaling the constant a, the constant 
To can be chosen equal to that in eq. (119). tq will be flxed below in terms of the initial condition 
z; the constant a on the other hand is flxed by matching with the gaussian solution. Notice that, 
as before, a does not depend on z. We do not know the explicit value of a, but we can argue that 
(J ~ 0{1). Indeed, the gaussian solution breaks down when t + Tq ~ 0{1), afterwards the solution 
will reach the linear regime in a time At ~ this means that the linear approximation breaks 

down when r + Tq ~ C^(l), implying a ~ 0{1). 

Let us now start the computation by flxing the relation between tq and z via 

= /ii„(0; z) = 1 - (76^° cos (Vero) z^ ae'° cos (Vero) , (121) 

which is vahd as long as Tq <S — 1. Since the solution / is constrained to be between and 1, in 
particular we need that for all r in < r < — Tq 

/iin(T;^)<l, (122) 

which gives the lower bound 

r„>-^. (123) 

Let us now study the analytic structure of /(r; 2;). Continuity in e suggests that /(r; 2;) has a 
branch cut at large enough negative z as before. As e approaches zero the origin of this cut Zcut 
goes to 2; = as well. To follow how Zcut depends on e at small e note that eq. (120) implies that 
/iin(T; z) is a function of z through its dependence on tq, upon which it depends analytically. So, 
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the only non-analyticity in z of /(r; z) can come from a non-analyticity of ro(-2). The boundary 
condition (121) tells us that z{tq) is analytic. In inverting this relationship, therefore, the only 
non-analyticity can arise if dz/drQ vanishes at some value of z, 

Q ^ o-e^° (cos(Vero) - Vesin(^To)) ^ ^ 

dTQ 1 - ae^o cos{^/eTo) 

=> cos{^/eTo) — ^/e sm{^/eTo) = 

^ -.--^-1. (125) 

By plugging this value into (121), we obtain 

z,^, ^ , (126) 

where we defined 

V; = . (127) 

With this in mind, we can begin to evaluate the probability distribution using the saddle point 
approximation. As in the previous section 

p(V,T) ^ ^ e-g(^o) = ^e-^(^°) , (128) 

S{z) = ^{r + Tof-zV, (129) 
where the saddle point Zq is determined by 

S'{zo)^l[r + ro{zo)]T;,{zo)-V^O ^ . ^ ■ (130) 

2 2a [cos( Ve^o) - ve sm{y/eTo)\ 

The above relationship implies that for volumes that are large but smaller than V = e"^, Tq is large 
and negative, so that the linear approximation in (121) is justified. t + Tq is also large and positive, 
so that we can trust the gaussian approximation for /(r; 2:0). As long as Tq > — tt/ {2^/e) = — log V^, 
the denominator does not vanish, and tq is approximately given by 

To^-logV. (131) 

The condition that tq > — log then reads V <Ve. When plugged back into eq. (128) the above 
solution gives the same form for p(V, r) as in the case fl > 1, namely 

p ^ 7\Ae-3(^-i°s^)' for K V < Min{F, 1/J , (132) 

where 
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There are two interesting cases: V and V > V^. Let us start with the case V < K- As V 
reaches V, tq -\- r becomes of order one, and we can not trust anymore the gaussian solution. For 
this reason, in exploring the regime V >V, analogously to what we did in the case of > 1, we 
close the contour along the cut in the negative real z axis (see fig. 11), and perform the integral 
of the imaginary part of /(r; z)e^^ 



f + OO 

-V\z\ 



p(V^,t) = ^ / ""^1^1 2ilm[f{T;-\z\)]e- 



(134) 



Because of the exponential suppression, the integral is dominated by \z\ ~ 1/V > |-2cut| for 
V < V^. In this regime, \z\ > |^cut|, and from eq. (121) we have tq ~ log{z/a). The condition 
\z\ ~ 1/y amounts to Re[ro] ~ —logV (notice that Im[ro] ~ iTi). Therefore in this regime 
T + Re[To] ~ \og{V /V) < —1, and we can use the linearized approximation /iin(T; 2;) for f{T;z). 
The imaginary part of fiin{r; z) on the negative z axis is rather complicated. However, we can 
approximate cos(v^(r + tq)) with a constant of order one. Then, the integral (134) can be 
estimated as 

p{V, t) I d\z\ e^z e-^l^l ^ , V <V <V, , (135) 



Zcut 



where we ignored constants of order one. We see that in this interval of volumes p{V, r) decreases 
as l/V^, exactly matching the analogous regime we found for > 1. 

As V keeps increasing and becomes larger than K, the solution to the boundary condition 
eq. (121), with z ~ 1/^, gives Tq ~ — 7r/(2y^) = —logVe- In this case, the integral on the 
discontinuity becomes dominated by the beginning of the cut, and we obtain 

piV,T) ~ / d\z\ e^e-^l^l ~ ^-e-^^V-eV/v. ^ V<V,<V, (136) 

where again we have ignored constants of order one, and where we have used that Zcut — 
— (Tv^/(eV^). For volumes larger than V^, p{V,t) decreases exponentially (this exponential tail 
was also recently found in [32]). 

Notice how the two solutions (135) and (136) glue together: for V < V^, Tq decreases like 
— logF giving a behavior to p(y); when the volume reaches ~ V^, tq freezes at a value 

~ (— logK) = ~7''/(2-\/e) and the exponential tail {e~^^^^^^) kicks in. Notice also that the 
point, where the exponential tail enters, goes to infinity for e — > smoothly matching the result 
for > 1. 

Let us now concentrate on the opposite regime. Namely, as e increases leaving r fixed, at 
some point becomes smaller that V. In this case, for V < V^, we still have the result of 
eq. (132). However, for larger V, tq does not continue to decrease as — logV, because in this case 
the denominator of eq. (130) goes to zero and determines the behavior of Tq. In particular we can 
expand the denominator of eq. (130) around Tq = — 7r/(2y^) 

2a [cos(y^ro) - V^sin(V^ro)] ^ la^e (^tq + ^ + 1^ , (137) 
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Figure 11: Contours of the anti-Laplace transform integrals in the z-plane (on the left) and probability 
distributions of the volume (on the right) for Q = 1 — e. 

First case \fe < tt/{2t) (first row); at small volumes the integral can be solved via saddle-point 
approximation and gives a gaussian-like distribution (in blue). Around V = V (gray region) our ap- 
proximations for the solution of the differential equation break down. At larger volumes the contour of 
integration can be deformed around the cut. As long as V '^V^ the integral is dominated by a region that 
is much larger than the distance Zcut between the beginning of the cut and the origin, the integral is thus 
equivalent to the integral over a cut that starts at z = (dotted red contour) giving a power-law behavior 
for p{V,t) (in light red). At V > V^, the contour integral "sees" the distance Zcut between the beginning 
of the cut and the origin, and it produces an exponential tail (in dark red) . 

Second case \fe > 7r/(2r) (second row).- this time the saddle point works for all the volumes; for 
V <Ve it produces a gaussian-like tail that is converted into an exponential tail for V >Ve. 

that indeed vanishes for tq = — 1. We can now substitute eq. (137) into eq. (130) and get 



where the third term in the expression for tq is small for V ^ V^. We see that Tq approaches 
asymptotically from the positive side the value tq = ~ ^- ^^is regime, r + tq is always 

larger than one and positive, and the gaussian approximation to our solution holds. Further, 
by substituting (138) in (130), it is straightforward to see that zq moves on the real axis and 
approaches ^cut from the positive side, reaching Zcut only as F ^ 00. This tells us that for V, < V, 
unlike in all the former cases, the saddle point approximation holds for all V^'s. Plugging in the 
expression for the saddle point integral we get 

p(V,r) « AAe"H^-27i)'-?^^^/^^ , V>Ve, (139) 

with 



3/2 /y\3/2 



N 3/2 
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We still have an exponential tail, with the same behavior as we found in the case in which V 
was smaller than V^. This is a very intuitive result: as e increases, the exponential tails becomes 
relevant at smaller and smaller volumes, and it eats away the main part of the distribution (see 
fig. 11). This also offers a consistency check between the two ways in which we are computing the 
inverse Laplace-transform: the saddle point approximation and the integral along the cut. 

Finally the normalization of p{V, r) changes smoothly when crossing the critical point 0,-1 

Pext = /(r; 0)^l-a^/~eTe'-^^ ^1- 6a^eN,-— . (141) 

In the limit e — > the normalization goes to one, but when the exponential tail starts to remove 
the bulk of p{V, t) (i.e. ~ = e^^'^) the volume normalization quickly drops to 0. 

For e — 0{l) our approximations break down and the calculation becomes more complicated: 
the bound on tq in eq. (123) becomes of 0{1) and both the saddle-point approximation and 
the integral along the cut arc dominated by the region of z where the linearized approximation 
no longer applies. However we can still say something about the normalization of the volume 
distribution Pext- For large r the solution will still be in the gaussian regime, and we can write 
f{T]z) as 

f{T;z)^k{Q,z)e-'-^^ , (142) 
where k and Tq are unknown constants that depend on ft and z. From eq. (46) we know that 



ext 



r p{V)dV = /(r; 0) = k{n, 0)6-^"^^^ . (143) 
Jo 



This formula does not tell us the explicit dependence on Q of Pext but since r we can 

extract the dependence of Pgxt on the classical number of e-foldings Nc (equivalent to the initial 
position of the inflaton) . For large Nc we get 

Pext ~ e-"(^^=)' = e-"('°^^=)' . (144) 
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The probability not to eternally inflate when il < 1 goes to exponentially with the square of the 
classical number of e-foldings. 

Finally, note that the exponential behavior of p{V) at V > implies that the moments of 
the distribution p{V) do not diverge at Q < 1. Naively, this disagree with the results of [16]. 
However, there is no conflict. The point is that here wc calculate moments using the probability 
distribution p{V) obtained by taking the infinite time limit. At < 1 its normalization is less than 
one, indicating that there is a contribution localized at infinity that is not taken into account (cf. 
with the discussion of the two-site example in section 2.1.1). On the other hand, in the calculation 
of [16] one first finds the moments and only then takes the infinite time limit. The latter procedure 
effectively takes into account the contribution at infinity producing diverging moments at Q < 1. 



3.5 $1 = 0: Deeply inside eternal inflation 

We are finally led to study the extreme hmit of eternal infiation: the case Q = 0, which corresponds 
to a completely flat inflationary potential. 

The job is quite easy in this case: the differential equation for / simplifies and we can find an 
explicit solution. With vanishing fl the solution to eq. (37) that stops on the top of the potential 
in an infinite time reads 

/(T;^) = e5-i(-+^™)'. (145) 

Then the probability distribution is 

p(V, r) = / dz e^i(-+^) , (146) 

and to evaluate this integral we can use the saddle point method as before (actually, since the 
integral is gaussian in V2 + 4:Z, this procedure is exact). 
The saddle point condition in this case reads 

S{z) = + i (r + ^2 + 4^) '-zV, (147) 



so that 



V V -1 

s"(«o) = (151) 



and the expression for the probability distribution is 



V-l 
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The tail of the distribution is again exponential, and this nicely matches the Q < 1 case we studied 
before. 

We can easily also calculate the normalization for p{V,T), 



which matches with the approximate formula of the previous section. As long as f2 ^ 1 the 
corrections from the friction term in the differential equation are small, so that eq. (152) is a good 
approximation in this regime. 

3.6 Realistic models: finite barrier effects and slow roll corrections 

So far, we worked in the approximation of an infinitely long inflaton potential and treated Q and H 
as constants. In a realistic situation both these assumptions do not hold: fl and H change slowly 
as functions of the inflaton field and the latter can vary only in a finite range. This may be, for 
example, a consequence of either quantum gravity effects, if the potential grows monotonically up 
to high vahics of the field, or a reheating region if the potential has a maximum, or the steepening 
of the potential itself at a certain region. Let us discuss at the qualitative level how these effects 
change our results. 

Let us begin with the consequences of the finite range of the infiaton field. The details of the 
underlying physical mechanism are not relevant for our qualitative discussion, so we introduce 
this effect by including a reflecting barrier in the stochastic process at large values of the inflaton 
fleld. In terms of the mechanical problem this implies that we are now looking for a solution that 
stops at a flnite time Tb (see (42)). 

The presence of a barrier affects our results in two different ways: it changes the tail of 
the probability distribution and it slightly shifts the critical value Qc for the transition to the 
eternal regime. Let us start with the first effect. In the non-eternal regime, we have seen that 
the probability distribution is peaked around a volume of order e"^^, where is given by (6) 
and changes between Nc for large Q and 2Nc at the transition point Q = 1. Given the relation 
Nc — T I (6\/n) this implies that the typical trajectory undergoes field excursions at most of 
order r up the infiaton potential. Hence putting a barrier at t\, does not affect the bulk of the 
trajectories, and therefore the probability distribution as long as the barrier is far enough from 
the starting point r <^ rj,. Still, the barrier cuts the trajectories that would have otherwise crossed 
it, and therefore we expect an additional suppression of the tail of the distribution for volumes 
corresponding to r > r^-, that is for N > N^-, where Ni, = Tb/ {Q^/Vt). 

As a result, for N > Nj, the probability distribution becomes exponentially suppressed as a 
function of the volume as opposed to the v~^~'^+^^- behavior that we found. There are several 
ways to see that the barrier indeed leads to the exponential suppression. For instance, as we saw 
in section 3.1, the function f{z) is analytic at 2; = when the barrier is present. Consequently, the 
whole integration contour for the Laplace transform can be deformed in the region with Ke{z) < 0, 
and at the large volumes the integral in (44) decays as e~^^(^="*)^, where Zcut is the singularity of 
the function / at the smallest distance from the real axis. 




(153) 
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In fact, recently an exponential tail of the volume distribution in the eternal regime was 
calculated in [32]. By extending the analysis of [32] into the non-eternal regime we verified that, 
as a result, the exponential suppression sets in for a number of e-foldings of order Nh. We do not 
provide details of this analysis here, as it would require an extensive introduction into methods of 
[32] and would take us too far away from the main line of our paper (instead, in the concluding 
section 4, we will present an intuitive argument explaining the origin of this exponential behavior). 
This also agrees with what we found in the two site models, where the probability distribution 
decreases exponentially at large volumes both in the eternal and the non-eternal regime. 

Notice that, since Q > 1 imphes Nc < Sds/^'^ (see eq. (4)), in the same regime we also have 
an upper bound on Ni,, which has to satisfy the same condition Ni, < This also means 

that the behavior of p{V) for volumes larger than e^"^^/"^ is always exponential in the volume, i.e. 
~ g-constv^ This property of the probability distribution will be further discussed in the next 
section in connection with the bound on the number of e-foldings. 

The second effect of the barrier is related to the first. The presence of the barrier cuts out 
some of the trajectories producing the largest volumes and thus may slightly delay the entrance 
in the eternal regime to f2c < 1- This shift can be calculated in the following way. Since f2c is 
defined as the value of Vt where Pcxt starts deviating from one, and since Pcxt = /(t; 0), we need 
to study the solution / for 2; — > 0. Consider some fixed < 1, and let us recall that in this case 
the oscillator is not anti-over-damped (eq. (37)). Around / = 1, we can linearize the potential 
and obtain an oscillating solution with a period 

independently of the initial velocity, as the oscillations arc harmonic (sec the dotted lines of fig. 7b). 
Consequently, in the linear regime it takes an amount of time equal to T/4 for the solution to 
come at rest. The returning force for the actual unharmonic potential (47) is smaller than for 
the corresponding harmonic potential, so it always takes longer than T/4 to come at rest (and 
can take arbitrarily long, as at < 1 there exists a solution that stops at the top of the hill 
in an infinite time). Consequently, if the barrier is close enough, Tf, < T/4, there is no solution 
with the appropriate boundary condition, /(r;,;0) = (except the trivial solution /(r;0) = 1), 
even at Q < 1 and inflation is not eternal. This argument implies that the critical value for Q is 
determined by 

Note that for Vt < i7c, the volume where the barrier effect sets in and changes the tail of the 
distribution is of order V = e^^'' = e'^/2\/i-^^c_ ^j^g^^ jg always exponentially larger then the 

volume Ve, where the exponential suppression found in section 3.4 sets in. Consequently, in the 
eternal regime the barrier affects only the tail of the distribution that is already suppressed as 

g— const y 

Let us now discuss how the dependence of Q and H on the infiaton field changes our results. 
Let us sec first when this dependence become important. Throughout this paper we have been 
interested in what happens close to the eternal regime in the limit Mp\ 3> H. In other words we 
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have taken the hmit H/Mp\ 0, while keeping Q fixed. This hmit imphes that we have been 
working in the extreme slow roll regime. Indeed, in this limit one has, 

H ^H^ 

However, this does not imply that one can completely neglect the field dependence of Q and H as 
we have another large parameter — the variation of the infiaton field or, equivalently, the number 
of e-foldings. For instance, by taking the variation of the Priedmann equation, one has 

Ml^HAH ~ V'A(t) 

and the condition that H can be treated as a constant reads 

So all our results apply if one takes the limit H/M-p\ — ^ while keeping ^2 fixed and small 
enough such that the inequality (156) holds. Note that (156) does not prevent us from considering 
arbitrarily large Nc if H/Mpi is taken to be sufficiently small. 

Nevertheless, one may wonder what happens to the shape of the volume probability distribution 
for longer infiaton trajectories, such that the variation of H (and Q) has to be taken into account. 
In this case the coefficients in front of the different terms in the mechanical equation acquire a 
dependence [35]. This clearly may affect the details of the shape of the probability distribution. 
However, as long as H <^ Mpi we still expect these effects to be small, with p{V) still sharply 
peaked around the average number of e-foldings N, which takes values between and 2Nc. One 
can check this statement by using the techniques of section 3.1 to calculate the average and the 
higher moments of the distribution. We can consider, for instance, the linear differential equation 
(57) for the average in the generic case, where Q and H can be functions of the infiaton field 
(or equivalently of r). In the variable r, defined now as 

T = 6 f VndN^, (157) 



the differential equation for the average is the same as in eq. (57) with a small correction (of 
order if^/M|j) to the anti- friction term, which now depends on r implicitly via fl. Since in the 
non-eternal regime the anti- friction term is minimized at Q = 1 , by setting n — Iwe would obtain 
a trajectory that has a faster velocity at each moment of time in order to stop at the same moment 
as required by the boundary condition (59). It follows that the average is always smaller than e'^, 
which, from cq. (3), is smaller than e^'^^^'^. This implies that the bound < Sds/Q still holds. 

In fact, we can also find an approximate expression for the average volume in the general case. 
Namely, for non-constant Jl, we can replace (61) by 

lim {V) = e^^-'^^ (158) 

T(,-+0O 
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which is a good approximate solution as soon as ^ fl. This gives for the average number of 
e-foldings 

37V- Juj-dT. (159) 

Analogous arguments can be applied for higher moments; as a result we see that the distribution 
remains peaked around the average value that takes value between N^. and 2Nc. 



4 Discussion 

To summarize, in this paper wc calculated explicitly the probability distribution p{N) for the vol- 
ume of the Universe after a period of slow- roll inflation (as before, by the number of e-foldings N we 
understand one third of the logarithm of the total volume produced during inflation, N — ^ log V). 
Our results cover both the eternal and the non-eternal regime. Let us start this concluding section 
by summarizing these results and then by explaining how all different kinds of behavior we found 
for p{N) have an intuitive physical explanation. 

In general, the function p{N) exhibits three qualitatively different regions. Namely, in the 
non-eternal regime > 1 it is peaked at N N, where the average number of e-foldings N is 
given by 

2N 

N = ' . (160) 



For N < N it has a gaussian form. 



1 + Vi - ^-^ 



p{N) oc e , (161) 



with a width a given by 



cr^ = , . (162) 



The behavior changes at > where the probability distribution becomes exponential in A^ (or, 
equivalently, power-law in the volume V), 

piN) oc e-^"^(^+V^) = y-^K^+V^) . (163) 

Finally, if a barrier preventing the inflaton to take arbitrary large values is present, this power-law 
tail becomes further suppressed at larger volumes and turns into an exponential in the volume 

p{N) oc e-^^'"^ = e-^^ . (164) 

This behavior sets in at A" ~ A^, where A^^ is the number of e-foldings on the classical inflaton 
trajectory from the barrier till reheating. 

What changes in the eternal regime, f2 < 1, is that the exponential behavior (164) sets in at 
A^ ~ 7r/6-\/l — ^ even if the barrier is absent. If (1 — fl) is not too small this happens at N < N 
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so that the gaussian regime (161) interpolates directly to the superexponential one (164) without 
intermediate exponential behavior (163). 

It is rather straightforward to understand the origin of the above three different types of 
behavior for p{N) at the intuitive level. First, to produce a small number of e-foldings N <^ Nc, 
the inflaton field during the first few e-foldings needs to perform a jump in the whole volume to a 
value of (p corresponding to an average number of e-folding equal to A^. The inflaton fluctuations 
away from the classical trajectory follow, at early times, a gaussian distribution (because of the 
random walk) and the size of the jump in fleld space is directly proportional to (A^ — A^), so the 
probability of such a jump is suppressed by e'''^^^~^^'^ (c is a constant) in agreement with our 
result (161) at small volumes. 

On the other hand, the least expensive way to produce a large number of e-foldings N ^ 
is for one Hubble patch to go high up the potential, till values of the fleld corresponding to have a 
classical trajectory with N e-foldings. If such a fluctuation happened, the probability to produce 
at least N e-foldings becomes of order one. To estimate the probability of such a fluctuation, 
note again that the probability distribution for the inflaton fluctuations A0 around the classical 
trajectory is gaussian. The crucial difference with the small volume case is that now it is not 
necessary for the fluctuation to happen in a short period of time, and the variance of the inflaton 
distribution grows linearly as a function of time t. As a result, the probability p of the fluctuation 
is maximized for times t corresponding to order N e-foldings, and depends exponentially on N , 
p oc e~'^i^'^ oc e~'^'^^ in agreement with (163). This argument is essentially identical to the one 
provided in [16] to explain why high enough moments of the volume distribution diverge in the 
non-eternal regime if no barrier is introduced. 

The last argument fails at sufficiently large volumes both in the presence of a barrier at high 
values of the inflaton fleld and in the eternal regime. In the flrst case it fails because there is 
a limit on the maximum length of the classical trajectory, while in the second because even if a 
fluctuation to high values of the inflaton fleld happened, one is not guaranteed to end up with a 
finite number of e-foldings (in fact, the higher the fluctuation is the smaller is the probability for 
inflation to terminate). In both cases one ends up with getting a superexponential suppression 
of the probability distribution at large volumes, (164). This can roughly be explained by the 
necessity for an exponential number of independent and rather improbable events (corresponding 
to the individual Hubble volumes produced during inflation) to happen. Namely, in the eternal 
regime, inflation should terminate for ~ e^^ Hubble patches and all their children (while in the 
non-eternal regime the probability for this to happen is equal to one). In the non-eternal case 
with a barrier, for N ^ N^, an order ~ e^^ Hubble patches have to five an anomalously long 
time, corresponding to a number of e-foldings much larger than Nj,. 

To summarize, we see that our results on the shape of the probability distribution of the 
reheating volume not only pass a number of consistency checks with results obtained by different 
methods, but also can be understood in a rather intuitive level. We believe making the explicit 
form of the volume probability available may help to understand the geometry of eternal inffation 
and offers also a natural "theoretical" observable that can be useful in further studies of eternal 
inflation. 

A further motivation for doing our calculation was to establish whether a quantum version of 
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the bound 

SN < cSds , (165) 

exists and, if so, whether this bound can be made sharp, i.e. whether there is a concrete vahie 
for the coefficient c in (165). Our results provide a positive answer to both questions. Namely, we 
find that the quantum version of the bound (165) can be formulated as follows: 

The probability for slow-roll inflation to produce a finite volume larger than e^''-^^'^, where Sds 
is de Sitter entropy at the end of the inflationary stage, is suppressed below the uncertainty due 
to non-perturbative quantum gravity effects. 

Indeed, our results imply that an inflaton trajectory with more than 2Nc e-foldings and such 
that infiation terminates globally in the entire space is exponentially unprobable, ~ e~^. Given 
that in the non-eternal regime we have the bound (4), we conclude that the probability for slow- roll 
inflation to last more than Sds/Q e-foldings and to terminate globally in the entire space is smaller 
than the uncertainty coming from non-perturbative quantum gravitational effects (~ e^'^'^^). In 
fact, there is an even stronger statement that follows from the observation of section 3.6 about 
the bound on the barrier position Nf,. Indeed, as long as ^2 > 1, A^^ is smaller than Sds/^'^i 
which implies that the transition to the regime where piV) drops exponentially with the volume 
^g-const-y-^ starts exactly when the produced volume starts violating the bound, V > e^'^^l'^ > e^^\ 
Consequently, the probability to produce much larger volumes V ^ e^ds/'^ jg super-exponentially 
suppressed, ~ . In turn this produces a super-exponentially small probability to produce 

volumes V ^ e^'^^^'^. 

We believe this provides a non-trivial test of the de Sitter complementarity idea. However, 
one may wish to go further and ask whether the particular value c = 1/2 that we found provides 
any insights into how the de Sitter complementarity works. Of course, all our calculations were 
done in the limit where gravity is non-dynamical, therefore it is not clear whether the value 
c = 1/2 really tells us something fundamental about the properties of de Sitter space. Let us 
however assume that it does and speculate on what would be the interpretation of this particular 
value. Recall that the original motivation for the bound (1) is coming from the idea that the 
black hole complementarity applies in the de Sitter case as well, so that the global effective field 
theory description of the FRW slices breaks down and information about the outside observers 
gets released in de Sitter fiuctuations. Then, if infiation ends, there is the danger of violating the 
linearity of quantum mechanics by creating a "quantum xerox" — one can see the information twice, 
first holographically in the de Sitter fiuctuations and later in a direct way when the corresponding 
mode comes in. 

Recall that the analogous paradox could potentially arise in the black hole case after one 
measures more than Shh/2 Hawking quanta {Sbh is the black hole entropy). In this case, the factor 
1/2 arises because, as pointed out by Page [38], if one measures k degrees of freedom of a larger 
system described by n degrees of freedom, generically, the resulting density matrix looks thermal 
and carries less than a single bit of information as long as A; < n/2. 

This does not give rise however to any paradox: indeed, if one waits a long enough time 
outside the black hole, so that ~ Shhf^ Hawking quanta are emitted, it is impossible to observe 
the same information a second time inside the black hole, since it gets destroyed by the curvature 
singularity. 
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It is tempting to interprete our result by saying that the de Sitter analogue of one Hawking 
quanta is produced every time a new Hubble patch is created from the volume where the observer 
lives {i.e every 1/3 of e-folding). Then the quantum-xerox problem would arise after Sas/G e- 
foldings but it does not because of the bound. 

Although it may sound quite natural, this interpretation raises a number of issues. First, as 
argued in [39], a single observer in eternal de Sitter can never run into the xerox paradox, as 
the largest amount of entropy that can be stored within a single causal patch is bounded by the 
Bekenstein area of the largest black hole that fits in a single patch. The latter is equal to Sds/^i 
and smaller than the minimal amount {Sds/'^) required to extract the first bit of information. 

However we are in a rather different situation here, as inflation eventually terminates. Therefore 
should the bound be not there, one would have the possibility to measure Sds/4: quanta, give them 
to a friend who leaves for a different Hubble patch, then measure other Sds/4: quanta and, after 
inflation ends, the two could meet and in this way collect Sds/'2 quanta in total. Then the bound 
suggests that for the purpose of extracting information it is legitimate to collect quanta in one 
Hubble patch and keep them until the end of inflation in another. However, this reasoning also 
implies that it should not be possible to use quanta collected in different Hubble patches for 
extracting information — otherwise one would be forced to conclude that the actual bound on the 
number of e-folding is much stronger, N < log S. 

Notice also that if inflation does not terminate everywhere, then there is no problem of duplica- 
tion of information if one observer lasts in the inflationary phase for longer than Sds / 6 e-foldings: 
when the observer undergoes reheating and enters the Minkowski phase she is never able to see 
all the rest of the de Sitter space. 

Another puzzle with drawing the parallel between the factor 1/2 in our bound and the one in 
the Page argument is that, if there are n hght species around then, hkely, the rate of how fast the 
information gets released is proportional to the number of species (at least this is so in the black 
hole case). On the other hand our bound does not depend on n (it could if the light species are 
scalars but does not for other fields). 

Finally, it is worth mentioning yet another reason why one should be cautious in taking the 
value c — 1/2 too seriously. Namely, here we focused on slow-roll models of inflation. In the 
effective fleld theory language of [37] these correspond to the models where the inflaton pertur- 
bations have sound velocity Cg equal to one. The classical version of the bound (1) was proven in 
[22] for general inflationary models, including those with small c^. In fact at small Cg the classical 
bound becomes even stronger, Nc < c^Sds- However the reason it arises is somewhat different. At 
small sound velocities the strong coupling regime, where the effective fleld theory breaks down, 
sets in always before the eternal inflation regime. So the bound (1) in this case follows from 
the requirement that the system is weakly coupled and the null energy condition is not violated. 
Consequently, at the present stage we cannot exclude the possibility of violating the bound in the 
strongly coupled regime, although this possibility appears highly unlikely. 

To summarize, we proved the quantum version of the bound on how long slow-roll inflation 
can last without becoming eternal. The existence of such a bound (eq. (165)) provides non- 
trivial support to de Sitter complementarity ideas, while it is still unclear whether the value of 
the coefficient c = 1/2 appearing in the bound that we found has a definite physical meaning. 
Before finishing, it is worth stressing that we believe that, independently of the answer to the last 
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question, explicit calculations providing a detailed quantitative understanding of the transition to 
the eternal inflation regime, like those we preformed in the current paper, have an independent 
value. Indeed, it appears to us that a detailed understanding of the geometry and the dynamics of 
the eternally inflating Universe might be an important step to reach a final verdict on the puzzling 
issues raised by the observation of the cosmic acceleration and the string landscape. 
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A Volume average from the inflaton stochastic equation 



As a cross-check of the method used in the paper, derived from the bacteria model, in this 
appendix we present an alternative calculation of the average volume (V) directly from the inflaton 
equations. The result is an improvement of a very similar computation performed in [16]. By 
definition the average of the volume is given by 



where the adimensional variable V is the volume expressed in units of the initial volume. This 
reduces to compute 



where Pr{t) is the probability that at a given point x the field reaches the reheating value 0,. at time 
t. By translational invariance this probability does not depend on the point x. Now the problem 
is reduced to compute Pr{t). The latter is related to the probability P{(j),t) for the infiaton to 
have a value at time t by 

Prit) = -jJ _d<pP{<p,t). (168) 
P{4>:t) can be found by solving the classic stochastic diffusion equation [34, 10, 11] 




(166) 




to 



(167) 



(169) 



where P{ip,a^) = P{4>,t) and 



= (p — (p — (pt 



(170) 
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is a Gaussian field with variance that grows Unearly with time 

.= = ^*. (171) 

■0 represents the fluctuations around the classical motion and undergoes a random walk. In the 
case the inflaton lives in a infinitely long potential and the reheating point is at 0^ — 0, the 
solution of (169) is given by 

P(^, a-) = (e-^ - e-^'^^Me-^) , (172) 

which implies 



Pr{t) = V2i^e-^'--^ . (173) 



We see that the integral in (167) converges only when 

^^-^|,>1. (174) 

This is the first signal that the system has entered the eternal inflation regime. For Q > 1, we can 
explicitly perform the integral in (167), to obtain: 

(V) ^ , ,.v..(vn-vn^)a ^[^^J^>>\ . n>l. (175) 



'67r2 



where 

is the volume in the classical hmit. The result, which agrees with eq. (118), is plotted in fig. 8, 
where we can see the singular behavior at Q = 1. 
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